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Abstract 


Part I of this report continues the investigation, initiated in pre- 
vious reports, of scattering from rectangular plates coated with lossy 
dielectrics. The hard polarization coefficients given in the last report 
are incorporated into a model, which includes second- and third-order 
diffractions, for the coated plate. Computed results from this model are 
examined and compared to measured data. A breakdown of the con- 
tribution of each of the higher-order terms to the total RCS is given. 
The effectiveness of the UTD model in accounting for the coating effect 
is investigated by examining a Physical Optics (PO) model which in- 
corporates the equivalent surface impedance approximation used in the 
UTD model. The PO, UTD, and experimental results are compared. 

Part II of this report presents a radar cross section (RCS) model, 
based on Physical Optics (PO) and the Method of Equivalent Currents 
(MEC), for a trihedral corner reflector. PO is used to account for the 
reflected fields, while MEC is used for the diffracted fields. Single, 
double, and triple reflections and first-order diffractions are included 
in the model. A detailed derivation of the Ej-polarization, monostatic 
RCS is included. Computed results axe compaxed with finite- difference 
time-domain (FDTD) results for validation. The PO/MEC model of 
this report compares very well with the FDTD model, and it is a much 
faster model in terms of computational speed. 
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I. HIGH-FREQUENCY TECHNIQUES 
FOR RCS PREDICTION 
OF PLATE GEOMETRIES 

A. INTRODUCTION 

Scattering from a rectangular plate and from simple targets composed of rect- 
angular plates are important and useful problems to study because specific scat- 
tering mechanisms can be isolated and closely analyzed. Some of the scattering 
mechanisms of current interest in developing high-frequency asymptotic modeling 
techniques are diffractions from electrically thin (thickness < A), lossy dielectrics 
backed by a perfect conductor; equivalent impedance representations for this ge- 
ometry; nonprincipal-plane scattering; and corner diffractions. These topics have 
been addressed in previous reports [1] - [8]. This part of the report details the 
continuing work on the problem of scattering from a rectangular plate backed by 
an electrically thin, lossy dielectric. 

To determine the effectiveness of the Uniform Theory of Diffraction (UTD) 
expressions for a coated wedge, these diffraction coefficients are used to model 
the principal-plane scattering from a strip/plate geometry coated with a lossy 
dielectric. The hard polarization case is examined in this report. Terms for higher- 
order diffractions are included in the model. A breakdown of the various terms is 
used to determine the contibution of the higher-order terms to the overall radar 

cross section (RCS). Computed results are compared to measured RCS data for 
validation. 
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One of the primary concerns in the modeling of coated conductors is the 
boundary condition describing the coating. This investigation has focused on 
using the Impedance Boundary Condition (IBC) and approximating the equiv- 
alent impedance of the coating backed by a perfect conductor in terms of a short- 
circuited transmission-line model [8]. This is considered to be an accurate approx- 
imation near and at normal incidence to the surface. In order to verify this, a 
physical optics (PO) model for the plate is compared to experimental data and to 
the UTD model. 

B. THEORY AND RESULTS 

1. Coated Plate UTD RCS Model — Hard Polarization 

In the last report [8], UTD coefficients for hard polarization scattering from the 
coated wedge of Fig. 1 were presented. These coefficients were derived from 
Griesser and Balanis’ UTD coefficients [9] for a general dielectric wedge. Adapting 
these coefficients for the special case of a coated wedge resulted in simplification 
of the expressions involved, which led to faster computational times. In a previous 
report [7], a model for the principal-plane RCS of the coated plate of Fig. 2 was 
presented. This model used only first-order diffraction terms. Results from this 
model, shown in the last report [8], revealed good agreement with experimental 
data near and at normal incidence to the plate. It was thought that the addition of 
terms to account for higher-order diffractions and surface waves would immensely 
improve the results away from normal incidence. 

Terms to account for second- and third-order diffractions have been added 
to the RCS model of the coated plate for the hard polarization case. The fields 
for the various orders of diffraction are given below. An e JU,t time convention is 
assumed and suppressed throughout. The incident field is: 
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Incident Field 


Ei = -a*£ 0 e i * (rco8 *' +! ' 8in *'> (1) 

The diffracted fields are formed in the usual way for a UTD solution as a product of 
the incident field, an amplitude spreading factor, a phase factor, and the diffraction 
coefficient. The necessary diffraction coefficients will not be repeated here as they 
can be found in the previous report [8]. As has been discussed in the last report, 
the equivalent impedance of the coated surface is modeled using a short-circuited 
transmission-line approximation. The appropriate equation is given also in the last 
report [8]. 

The first-order diffracted field is: 

First-Order Field 


D-jkP 


Edi = -a+E 0 Vk-j=r [e-^^' +c ^D pwff (0;,0 1 X,2) 




\f£~P 

,+co ^ ) T> w/ ( 0 ', ^2X2)] 


( 2 ) 


^pw// o> n ) refers to the diffraction coefficient for plane- wave incidence, 

far-field observation. The explicit expression is given in Eq. (15) of the previous 
report [8]. This coefficient is analogous in form to the Keller GTD diffraction 
coefficient for perfectly conducting wedges. The parameters needed in Eq. (1) are 
defined in Figs. 1 and 2, where Oq is the Brewster angle for the coated side of the 
plate. This is described in Section II. A of [8]. 

The second- and third-order diffraction terms are expressed in terms of 
Dcwff and D vru ji (d, 0', 0,#o> n )> the diffraction coefficients for 
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cylindrical- wave incidence, far-field observation and for plane-wave incidence, ob- 
servation at a finite distance, respectively. The expression for D cwJf () results from 
the appropriate modification, discussed in Section II.B.3 of [8], of Eq. (16) of [8]. 
Dp W fd ( ) is given by Eq. (16) of [8]. The second- and third-order fields are: 

Second-Order Field 


E d2 = -a '. < j > E 0 '/k 


e -jkp e -jkw ^ 


x e J 2 


y/ltp y/w 

(m,V>;,0°,9‘,2) D„„ (w, 0°,*,#‘,2) 


+ 2) 0”,*,«},2)j (3) 

Third-Order Field 


Ed3 


—a 'jEoVk 


e~ jkp e -ju™y/k 
y/Vp W 


x Dp W /d(w ,0 ,0 , > 2 ) Dcwfj (w, 0°, 0°, 6q , 2^ 

X («,, 0°, 01, 0$, 2) 


+ 


e i“ i (cOS ^'+COS 


^Dcwff (ty,O o ,02,^> 2 )] 


(4) 


The RCS, computed using the preceding fields, is shown in Fig. 3 for a square, 
coated plate with a width of 2A. The coating has thickness t = 0.0423A with ma- 
terial parameters /i r = 1.539 -jl.2241 and e T = 11.826 — j'0.16639. The frequency 
of operation is 10 GHz. Comparisons axe made in Fig. 3 among the measured 
RCS, the RCS computed using the first-order field of Eq. (1) only, and the RCS 
computed using the fields of Eqs. (l)-(3). It is obvious that the addition of higher- 
order diffractions did not improve, or even noticeably change, the modeled RCS 
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values. The reason for this is shown in Fig. 4, which shows a breakdown of the 
field contributions. The second- and third-order fields are much smaller than the 
first-order fields; therefore, they are negligible. During the next reporting period, 
the surface-wave terms will be added to the model to see if these terms improve 
the model. 


2. Coated Plate PO RCS Model 

Thus far, the analysis of the coated plate incorporates a model of the impedance 
of an electrically thin coating backed by a perfect conductor that uses a short- 
circuited transmission-line approximation. This is considered a valid and conve- 
nient approximation for this geometry [10]; however, the peak experimental and 
theoretical UTD RCS values on the coated side of the plate do not agree exactly 
for the example considered in Fig. 3. Reasons for this could be the difficulty in 
accurately determining the constituitive parameters of the coating material and 
the extreme sensitivity of the theoretical model to even very small changes in the 
values of these parameters. Another reason could be inaccuracy in the UTD model. 
In order to validate the UTD model at this point, a PO model that incorporates 
the short-circuited transmission-line approximation for the coating impedance was 
derived for the plate; and the results were compared to the UTD and experimental 
results. 

The PO model was derived in the usual manner [11]. The resulting monos- 
tatic RCS equations for the perfectly conducting and coated side of the plate are: 


Perfectly Conducting Side 


<7 = 4t 


\ A / l kw sin 6 ) 


( 5 ) 
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Coated Side 


<7 = 4k 



sin(fcu;sin0)\ 2 
kw sin 9 ) 


( 6 ) 


where w is the plate width, L is the plate length, k is the propagation constant 
in free space, and 9 is the incidence/observation angle measured from the nor- 


mal to the plate. T r is the reflection coefficient [11] at a material interface with 
an impedance, Tj eq . For the results in Fig. 5, Tj eq is taken as the short-circuited 

transmission-line equivalence. The expression for this is given in Eq. (1) of the last 
report [8]. 


The results in Fig. 5 indicate that the UTD solution incorporating the short- 
circuited transmission-line approximation is accurate at normal incidence because 
the model agrees fairly well with the PO model at normal incidence to the coated 
side of the plate (<f> = 90°). The two models differ by approximately 0.84 dB 
at this point; however, the experimental results at this point are 1.15 dB higher 
than the PO results, indicating that most of the difference between the UTD and 
experimental results is due to inaccuracy in characterizing the coating material. It 
is encouraging to note that the UTD model agrees with the experimental results 
much better than the PO model in most areas. 


C. FUTURE WORK 

Future work will focus on three areas — the coated plate, nonprincipal-plane scat- 
tering, and the coated dihedral corner reflector. The analysis of the coated plate 
will be completed so that the UTD model includes both polarizations, surface-wave 
terms, and surface-wave transition region terms. It may be necessary to explore 
the use of higher-order boundary conditions to obtain accurate results at scattering 
angles that are not in the vicinity of normal incidence to the plate. Another area 
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of consideration will be accurately predicting the scattering at and near grazing 
incidence, a region which involves overlapping transition regions in which the tra- 
ditional application of UTD fails. 

Nonprincipal-plane scattering from both perfectly conducting and coated 
plates will be investigated using the Method of Equivalent Currents (MEC) and 
hybrid techniques. Much work has already been presented on nonprincipal-plane 
scattering from perfectly conducting plates in previous reports. Future work will 
involve completing the perfectly conducting plate model by incorporating corner 
scattering, either in terms of a revised equivalent currents model or through the 
use of a hybrid technique such as a combination of MEC and Moment Method 
(MM) or of the UTD and the Finite-Difference Time-Domain (FDTD) technique. 
Nonprincipal-plane scattering from a coated plate will be investigated using newly 
developed Incremental Length Diffraction Coefficients (ILDC) for the dielectric 
wedge [12] and other techniques. 

The primary objective of investigating plate scattering is to refine modeling 
techniques for scattering configurations that are part of more complicated target 
geometries so that the scattering from complex targets can be more completely 
understood and predicted using high-frequency techniques. To this end, the work 
that has been done on scattering from the coated plate will be incorporated into 
modeling the RCS of a coated dihedral corner reflector. This geometry will even- 
tually serve as a building block in even more complicated geometries. 
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Figure 1: Impedance wedge geometry. 


Plane-wave 

Incidence 



Figure 2: Geometry for principal-plane scattering from a strip/plate with a 
finite-thickness coating backed by a perfect conductor. 
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Figure 3: Principal-plane monostatic RCS of a coated plate (w = L = 2A, / = 
10GHz, coating: t = 0.0423A, fi r = 1.539 - j 1.2241, e r = 11.826 - jO. 16639). 



Figure 4: UTD components of the monostatic RCS of a coated plate (w — L = 2A, 
/ = 10GHz, coating: t = 0.0423A, (i r = 1.539 - j 1.2241, e r = 11.826 - 7*0.16639). 
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Figure 5: Principal-plane monostatic RCS of a coated plate (w = L = 2A, / 
10GHz, coating: t = 0.0423A, fi r = 1.539 — j'1.2241, e r = 11.826 - jO. 16639). 
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II. A PHYSICAL OPTICS/ 

EQUIVALENT CURRENTS MODEL 
FOR THE RCS OF A 
TRIHEDRAL CORNER REFLECTOR 

A. INTRODUCTION 

Corner reflectors are very interesting radar targets because they provide a large 
bistatic or monostatic radar cross section over a broad range of observation angles 
and frequencies. The two most well-known corner reflectors are the dihedral and 
the trihedral corner reflectors. The dihedral is formed by the intersection of two 
flat plates, whereas the trihedral is formed by the intersection of three flat plates. 
The fact that the trihedral has one plate more than the dihedral is also the reason 
why its backscatter cross section is much larger than that of the dihedral. In 
addition, the trihedral has a large radar cross section for any incident angle 
whereas the dihedral has a large radar cross section only for the case where the 
direction of the incident plane wave is normal to the dihedral axis (#, = 90°). 

Several researchers in the past have analyzed the scattering properties of 
dihedral corner reflectors using either Geometrical Optics (GO) together with Ge- 
ometrical Theory of Diffraction (GTD), or Physical Optics (PO) together with 
Physical Theory of Diffraction (PTD) [13]-[15]. In all those cases, the direction of 
incidence was always normal to the dihedral axis and therefore the analysis was 
carried out as if the object were two-dimensional. On the other hand, the trihe- 
dral corner reflector requires a full three-dimensional analysis whose formulation 
is much more complex than that of the dihedral. 
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In this analysis, Physical Optics (PO) and Method of Equivalent Currents 
(MEC) were applied on a square trihedral corner reflector to evaluate its monostatic 
radar cross section. Physical Optics was used for the calculation of first-, second-, 
and third-order reflections from the trihedral plates, whereas Method of Equivalent 
Currents was used for the calculation of the first-order diffractions from the edges. 
It is important to mention that for a relatively large trihedral corner reflector 
the first-order diffractions are very small compared to the backscatter fields due 
to internal reflections; therefore, the diffracted fields can be considered negligible. 
However, as the size of the trihedral becomes smaller, the first order diffractions can 
contribute significantly to the total backscatter fields. Furthermore, in the analysis 
of the trihedral corner reflector, it was assumed that higher order diffractions, as 
well as diffraction-reflections and reflection-diffractions, are negligible. 

B. ANALYSIS AND RESULTS 

The backscatter cross section of the square trihedral corner reflector, shown in Fig. 
6, is evaluated by considering single, double, and triple reflections as well as first- 
order diffractions. The Physical Optics approximation was used for the calculation 
of the reflected fields in the far-field region. For the evaluation of the diffracted 
fields, the Method of Equivalent Currents (MEC), introduced by Michaeli, was 
applied on the trihedral edges. The entire analysis was carried out for the case of 
Eg polarization only. 

1. Physical Optics (PO) 

Single, double, and triple reflections from the perfectly conducting plates of the 
square trihedral corner reflector were determined using the Physical Optics ap- 
proximation 

Jpo = 2n x H‘ (7) 
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where n is the unit vector normal to the plate under consideration, and H‘ is the 
incident magnetic field. It is important to point out that the above PO model 
becomes more accurate as the dimensions of the object increase. After evaluating 
the surface current density on a particular plate, the vector potential is calculated 

by 


A 

L, 


in J Is Jp ° R dS ~ L r J J s 3 Po eJkL ‘d* 

x sin 9 a cos <j) a + y sin 9 a sin <f> a -f- z cos 9 a 


-jkr 


( 8 ) 

( 9 ) 


In the case of monostatic RCS, 9, = 9i and <f> a — <£,•. The above integral should 
be evaluated over the illuminated area of the plate. For single reflections, the 
illuminated area is the entire square plate, whereas for double and triple reflections, 
the illuminated area has the shape of a polygon that changes according to the 
direction of incidence. From the vector potential expression, the corresponding 
far-field spherical components of the electric field can be calculated using 


Eg = —juiAg and Ej, = —juAj, (10) 

The E r component is considered to be very small and is excluded from the calcula- 
tions. In most geometries that are used in scattering problems, it is usually more 
convenient to modify the above general expressions to the ones given below [11]. 


E T 

~ 0 

(11) 

Eg 

—jkr)Nge~i kT 
4tt r 

(12) 

E* 

—jkr)N^e~^ kr 
4tt r 

(13) 

Ng 

= J cos cos 4>a + Jy cos sin 4>s ~ Jz sin 9 a )e^ kL ‘ds 

(14) 

Nt 

= J f r S ’ D 4>s + Jy COS (f> a )e jkL, ds 

(15) 
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The actual expressions for Ng and N+ are different for each trihedral scattering 
component. These components, in the case of Eg polarization, include the follow- 
ing: 

Component 1 : Single reflection from plate 1 
The incident magnetic field on plate 1 is given by 


= (-«* sin fa + k y cos fa)H Q e jkL ' (16) 

Li = x sin 0{ cos <f>i + y sin 0, sin <f>i + z cos 0, (17) 

The current density on plate 1 due to the incident magnetic field is 

Ji = (—a* cos <f>i - k y sin fa)2H o e jk{xaia0 ' “■*+»*»* .in*,) ( lg ) 

The corresponding expressions for Ng and are 

Ng = —2H 0 abcos 0 3 (cos fa cos <f> 3 + sin fa sin <f> 3 )e^ x+y ^ (19) 

.sin X sin Y . 

x ( ^r )( — > 

N<j, = 2H 0 ab{s\n fa cos fa — cos <f>„ sin fa)e^ x+Y ^ (20) 

.sin X w sinT . 


( OJAX . / Olll 1 

— X— ) 

ka 

X = — (sin 9 a cos <f>, + sin 0, cos fa) (21) 

kb 

T = — (sin 9, sin <f> a + sin $i sin fa) (22) 

The integrals in equations (14) and (15) are evaluated over the entire plate, since 
it is totally illuminated. 

Component 2 : Single reflection from plate 2 

H* 2 = (-a* sin fa + a v cos fa)H 0 e }kLi (23) 
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J 2 = k,2H 0 cos(f> i e jk{ysinS ‘ sin<t, ' +zcoee ''> 


(24) 


N e = -2 H 0 bc cos fa sin S ' m ^ ) (25) 

1 Zj 

^ = 0 (26) 

The integration for the evaluation of N$ was taken over the entire plate. 
Component 3 : Single reflection from plate 3 

H3 = (— a z sin <f>i + k y cos <j>i)H 0 e jkLi (27) 


J 3 = k z 2H 0 sin ^e ? *( xsin9 ' coe ^+ zcosff i) ( 2 8) 

Ng = -2H 0 acsin<f>i sin^e j( ' Y+z) ( S ^^ )(— - ) ( 29 ) 

A Z 

N * = 0 (30) 

kc 

Z = — (cos^, + cos 0 {) (31) 


Like the two previous scattering components, the integration in equation (14) was 
evaluated on the entire plate in order to obtain Ng. 

Component 4 : Double reflection - Plate 1 to Plate 2 

The incident magnetic field on plate 1 is given by (16). Then, Geometrical Optics 
(GO) is used to find the reflected magnetic field from plate 1, which is now the 
incident magnetic field on plate 2. The use of Geometrical Optics for the first re- 
flection will simplify the final expressions for the doubly reflected fields. If Physical 
Optics were used to account for the first reflection, instead of Geometrical Optics, 
then the expression for the scattered field in the far-field region would include a 
quadruple integral with non-constant limits, which is much more difficult to eval- 
uate than a double integral. The incident magnetic field on plate 2, which was 
initially reflected from plate 1, is given by 

Hi 2 = (—a* sin fc + a,, cos (f>i)H 0 e lkLl * (32) 
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Z/i 2 — x sin 9, cos (f>{ -j- y sin $,■ sin 4>i — z cos $,• (33) 

The current density on plate 2 due to the above magnetic field is 


J 12 = a z 2H 0 cos6 t e jk< ' ys ' me ‘ sin *•' ~ z cos * ) 


(34) 


The above expression is valid only on the illuminated area of plate 2. The shape 
of the illuminated area is a polygon, shown in Fig. 7, whose coordinates can be 
found accurately, at any angle of incidence, using trigonometric identities. In the 
case of double reflection from plate 1 to plate 2 


N * = (35) 

The expression for Ng can be found after evaluating the following integral over the 
illuminated area shown in Fig. 7 (The shape of the illuminated area changes with 
respect to the incident angle). 


Ng = 

— 2H 0 sin <f>i sin 0 3 J J 

(36) 

W yl2 = 

sin 0 a sin (j> 3 + sin 6 , sin <f>i 

(37) 

w zl2 = 

cos 0 3 — cos 0, 

(38) 

Similar discussion applies for the remaining doubly reflected fields. To save time 
and space, we present only the necessary equations for those components. 
Component 5: Double reflection - Plate 1 to Plate 3 

Hi 3 = 

(—a* sin + a„ cos <f>i)H 0 e jkL 13 

(39) 

L\3 ~ 

x sin 0 t cos 4>i + y sin 0, sin <f>i — z cos 0, 

(40) 

3 13 = 

k z 2H 0 sin 4 c^^-zco •«.-) 

(41) 

Ng = 

—2H 0 sin 4>i sin 0, J J e^ k ^ xWxl3+zW ‘ 13 ^ds 

(42) 
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0 


(43) 


iV* = 

Wxi 3 = sin 9, cos <f>, + sin 0, cos fa 
W^i 3 = cos 9 a — cos 9i 

Component 6 : Double reflection - Plate 2 to Plate 1 
H 2i = (a* sin fa 4- a y cos fa)H 0 e jkIj2i 
L 21 = —x sin 9, cos fa + y sin 9 % sin fa + z cos $i 

J 21 = (— a x cos + a y sin C08 ^' +y sulS| s,n '^ , 4 

Ne = 2 /f 0 cos ^ S (sin sin — cosfacos<f> a ) 

x J J e i k (^x2i+y w V 3i)d s 

Nj, = 2 i/ 0 (cos <f>i sin <f> a — sin fa cos <f> a ) 

X J j ^W^yW yil ) ds 

Wxii = sin 9, cos <f> a — sin 8 t cos fa 
Wy 2 i = sin 9, sin <f> a + sin 9{ sin fa 
Component 7 : Double reflection - Plate 2 to Plate 3 

H 23 = (a* sin fa + a y cos fa)H 0 e jkIj23 


L 23 

J23 

N e 

N* 


= —1 sin 9i cos fa + y sin 0,- sin fa + z cos 0, 
— _£ 2 sin <f>i2H 0 e^ k ^~ xamSi cos ^»+* «*<*;) 

= 2 H 0 sin 9 a sin fa n e jk(zW x 33 +zW l 33 )d s 


= 0 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 
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W 9 23 

= sin 9, cos <j) a — sin 0, cos </>,• 

(58) 


w a 23 

= cos 0, + COS $i 

(59) 

Component 8: Double reflection - Plate 3 to Plate 1 


h* 31 

= 

(— a r sin fa — a v cos <f>i)H 0 e^ kL31 

(60) 

L 31 

= 

x sin 6i cos <f>i — y sin 9i sin <^», + z cos 9{ 

(61) 

J 31 

= 

(a r cos fa - k y sin 

(62) 

N e 

= 

—2 H 0 cos 9 a { sin fa sin <f> a — cos <f>i cos <f > a ) 

(63) 


X 

j J e i fc ( lW *3i +y^«2i ) ds 


N* 

= 

— 2H 0 (cos <j>i sin <j>„ -f sin <f>i cos <f> a ) 

(64) 


X 

j J e jk(xW l31 +yW ail ) ds 


w x 31 

= 

sin 9 a cos <f> a + sin 9i cos <f> t 

(65) 

^*31 


sin 9 a sin <f>, — sin 0, sin 

(66) 

Component 9: Double reflection - Plate 3 to Plate 2 



HU 

= (— a r sin <f>i — k y cos (f> l )H 0 e jkL32 

(67) 


^32 

= x sin 9i cos <j>i — y sin 9 , sin <f>i + z cos 9i 

(68) 


J32 

= -k z cos ^.2// 0 e j * ( - y8in * C08tfl) 

(69) 


No 

= 2 H 0 sin 9 a cos & J J e^ w ^ +zW ^ds 

(70) 


N* 

= 0 

(71) 


Wy32 

= sin 0 S sin <f> 9 — sin 6{ sin <f>i 

(72) 
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cos 9 , -f cos 9 { 


(73) 


W, 3I = 

The remaining scattering components of the square trihedral corner reflector, 
besides the edge diffractions, are the triple reflected fields. The triple reflected 
fields are very important scattering components because they contribute the most 
to the monostatic radar cross section. The reason is that after three consecutive 
reflections the reflected fields return directly back to the source. This property of 
the trihedral corner reflector can be clearly seen if someone applies Geometrical 
Optics on the three intersected plates and finds the direction of propagation for the 
triple reflected field. In this analysis, Geometrical Optics was implemented for the 
calculation of the triple reflections. More precisely, Geometrical Optics was used 
for the first and second reflections in order to find the expression for the incident 
field on the third plate; then Physical Optics was applied on the illuminated area 
of the third plate. As in the case of double reflection, the shape of the illuminated 
area of the third plate, see Fig. 7, can be determined accurately using trigonometry 
and vector analysis. Geometrical Optics approximation gives very good results for 
the monostatic radar cross section of the trihedral, but not very good results for 
the bistatic radar cross section because the sidelobes are not predicted very well. 
To get more accurate results for the case of bistatic radar cross section, Physical 
Optics should be applied for all three consecutive reflections, which is certainly 
not an easy task to carry out. The reason is because strict application of Physical 
Optics would require the evaluation of six-fold integrals with non-constant limits. 
The triply reflected fields of the square trihedral corner reflector are the following. 
Component 10 : Triple Reflection from Plate 1 to Plate 2 to Plate 3 
After applying Geometrical Optics for the first two reflections and the required 
boundary conditions for the electric field on the conducting plates, the incident 
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magnetic field on the third plate can be found. 


h 123 = (a r sin </>,• + a y cos <f>i)H a e jkLi23 (74) 

L 123 = —x sin 0i cos <f>i + y sin 0i sin <f>{ — z cos 0i (75) 

According to Physical Optics approximation, the surface current density on the 


illuminated area of the plate 3 is given by 

J 123 = -a,sin^2 J ff o e J '* ( - xsin<?>COB0 *- ZCO8<?| ) (76) 

The expressions for Ng and are the following. 

Ng = 2H 0 sin<f>sin 0 a J j e jk ( xWll23+zW * 123 Us (77) 

N+ = 0 (78) 

Wri 23 = sin 0 1 cos <f> 3 — sin 0{ cos <f>{ (79) 

W z \2Z = COS 0 t - COS 0i (80) 

Component 11 : Triple Reflection from Plate 1 to Plate 3 to Plate 2 

H i32 = (~a* sin fa - ky cos fa)H 0 e ]kLl32 (81) 

T 132 = x sin 0{ cos <f>{ — y sin 0{ sin <f>i — z cos 0{ (82) 

J 132 = —a* cos < f >i 2H o e^ k ^~ yain0iaia<t>i ~ zcoee ^ (83) 

Ng = 2H 0 cos fa sin 0* J J s e jk ^ vWyl32+zWz132 ^ ds (84) 

N+ = 0 (85) 

!T c i 32 = sin 0, sin <f>, — sin 0, sin fa (86) 

W z i 32 = cos 0, — cos 0, (87) 
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Component 12 : Triple Reflection from Plate 2 to Plate 1 to Plate 3 



= (a* sin fa + a y cos fa)H a e ikL213 

(88) 

L 213 

= —x sin 0, cos fa + y sin 0, sin fa — z cos 0, 

(89) 

J213 

= — a z sin fa2H 0 e J ^ ~ x sin 6 ' co * ^ _z co * e > ) 

(90) 

N e 

= 2H 0 sin fa sin 0, J J e +zW z2 i 3 ) 

(91) 

N* 

= 0 

(92) 

W x 2 i3 

= sin 0, cos <j> 3 — sin 0 , cos fa 

(93) 

W M 213 

= cos 0, — cos 0, 

(94) 


Component 13: Triple Reflection from Plate 2 to Plate 3 to Plate 1 


^231 

= 

(a* sin fa — a v cos <f>i)H 0 e jkL ’ 23 '- 

(95) 

L 231 

= 

—x sin 0i cos <f>i — y sin 0{ sin fa + z cos 

(96) 

J231 

= 

(a, cos fa + k y sin fa)2H 0 e jk{ ~ xaia9iCoa4,i - yain$ ' ain ^ 

(97) 

N e 

= 

2 H 0 cos ^(sin fa sin <f>, + cos fa cos <f> 3 ) 

(98) 


X 

J J e M*w xi31 +z\v y231 ) ds 


N* 

= 

2 i/ 0 (sin <(>„ cos fa — cos <f> 3 sin fa) 

(99) 


X 

J J e i k ( xW *™+zW v231 ) 


W x2 31 

= 

sin 0 S cos <j) 3 — sin 0 , cos fa 

(100) 

W z231 

= 

sin 0 3 sin <f> 3 — sin 0 , sin <f>i 

(101) 
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Component 14 : Triple Reflection from Plate 3 to Plate 1 to Plate 2 


H|n 2 

= (— a x sin (f>i — k y cos <j>i)H 0 e ikL 3,2 

(102) 

L 312 

= x sin 9i cos <j>{ — y sin 0, sin <f>i — z cos 0 , 

(103) 

J312 

= —k z cos <j>i 2 H 0 e ik( ~ ysin9isin,t>i - zcose ^ 

(104) 

Ne 

= 2 H 0 cos <f>{ sin 0 , J J e ^ k ^ yWy312+zWz312 ^ ds 

(105) 

N 4 

= 0 

(106) 

W y312 

= sin 9 a sin (j), — sin 0,- sin 

(107) 

w * l2 

= COS 9 a — COS 0, 

(108) 


Component 15: Triple Reflection from Plate 3 to Plate 2 to Plate 1 


H 32 I 

— 

(k x sin <f>i — a„ cos <f>i)H 0 e jkL321 

(109) 

L 321 

= 

—x sin 0, cos <f>i — y sin 0,- sin + z cos 0,- 

(110) 

J321 

= 

(a* cos + k y sin <^)2 H 0 e jk{ ~ xaineiCoa,t> '- ysiD8i sin *-) 

(111) 

N e 

= 

2H 0 cos 0 a (sin 4>i sin 4> a + cos cos <f> a ) 

(112) 


X 

J J e i*(*w,3j 1 +zW r v32 j)^ 5 


N+ 

— 

2// a (sin cos ^ — cos <f> a sin <f>{) 

(113) 


X 

J j e jk(xw xi21+z w t321 ) ds 


W x 32i 

= 

sin 6, cos <f> a — sin 0, cos <f> a 

(114) 

W z 32i 

= 

sin 9 a sin <f> a — sin 9i sin <f>i 

(115) 
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2. Method of Equivalent Currents (MEC) 

The Method of Equivalent Currents (MEC), introduced by Michaeli, was applied 
on the edges of the square trihedral corner reflector to find the first-order diffrac- 
tions. The expressions for the electric and magnetic currents at the trihedral edges 
are similar to those published by Michaeli [16, 17]. For completeness, these equa- 
tions are given below for the specific case of an edge parallel to the z-axis. Similar 
expressions can be obtained for the edges at other orientations. 

M S = H i [ sin ^ sm^TT-o^/iV) 

zo NksmPsmP' [sinai cos((7r - ai)/N) - cos(<t>'/N) 

sin(A r 7r - <ft) sm((n - a 2 )/N) 

sina 2 cos((7r - a 2 )/N) -f cos (<t>’/N) 
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(122) 


<*i = -j In(/ii + \fnl - 1) 

a 2 = -j \n(fi 2 + \Jnl - i) (123) 

The above formulations were applied on the trihedral edges. However, the diffracted 
fields are relatively small compared to the reflected fields; therefore, they will not 
contribute very much to the monostatic radar cross section. On the other hand, it 
is important that these diffractions be included for the bistatic radar cross section, 
since now the reflected fields will be of lower magnitude and first-order diffractions 
will be more significant. In our analysis, first-order diffractions were included for 
both monostatic and bistatic radar cross section, because first-order diffractions 
become increasingly important as the target becomes smaller; even for the mono- 
static case. 

3. Results 

A complex program based on the above formulations has been developed to calcu- 
late both monostatic and bistatic RCS of the square trihedral corner reflector for 
the E$ polarization. This program is applicable for any angle of incidence and/or 
observation angle. It also runs very well and fast on various computers such as the 
SUN, VAX, CMS, or even PC. Two plots for the monostatic case are presented in 
this report. Fig. 8 shows the monostatic RCS for an incident angle of <f> = 45° as 
the angle 0 varies from 0° to 90°. Fig. 9 shows the monostatic RCS for an incident 
angle 6 = 66° as the angle <f> varies from 0° to 90°. For both plots, the size of each 
square trihedral plate was taken to be 5A. FDTD data were provided to compare 
our results. Comparison between the two methods shows a very good agreement. 
In addition, it is important to note that running an FDTD program to calculate 
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the monostatic RCS of a trihedral definitely requires a lot of time, whereas running 
the “PO + MEC” program to do the same task takes only a few seconds. 

C. FUTURE WORK 


In this report, only the monostatic RCS of the square trihedral corner reflector in 
the Eg polarization was discussed. Actually, the formulations given in the above 
sections can be used for bistatic RCS without any modification. However, we do not 
present any plots for the bistatic RCS because of lack of accuracy at the sidelobes. 
There is very good agreement with FDTD data at and near the mainlobe but 
not very good agreement at the sidelobes. The reason is because the Geometrical 
Optics approximation was used for both double and triple reflections. To get better 
graphs for the bistatic case, it is necessary that we use strictly Physical Optics for 
double and trible reflections. However, this is a very challenging task, since it would 
be necessary to solve a quadruple and a six-fold integral with non-constant limits. 
This will be a priority goal for the future. In addition, the formulation for the E 0 
polarization will be investigated to have a complete picture of the backscattering 
of the square trihedral corner reflector. Lastly, a similar procedure will be applied 
for a triangular trihedral corner reflector. Measurements will also be performed to 
compare with predictions. 
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Figure 8: Monostatic RCS of the square trihedral corner reflector with dimensions 

a = b = c = 5.0A, incident angle <f> — 45°, and polarization of the E-field in the 
theta direction. 



Figure 9: Monostatic RCS of the square trihedral corner reflector with dimensions 
a = b = c = 5.0A, incident angle 9 = 66°, and polarization of the E-field in the 
theta direction. 
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